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The classical region of a Bose gas consists of all single-particle modes that have a high average occupation 
and are well-described by a classical field. Highly-occupied modes only occur in massive Bose gases at ultra- 
cold temperatures, in contrast to the photon case where there are highly-occupied modes at all temperatures. For 
the Bose gas the number of these modes is dependent on the temperature, the total number of particles and their 
interaction strength. In this paper we characterize the classical region of a harmonically trapped Bose gas over 
a wide parameter regime. We use a Hartree-Fock approach to account for the effects of interactions, which we 
observe to significantly change the classical region as compared to the idealized case. We compare our results 
to full classical field calculations and show that the Hartree-Fock approach provides a qualitatively accurate 
description of classical region for the interacting gas. 
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I. INTRODUCTION 

There has been much recent work on using the classical 
field approximation to model Bose-Einstein condensates at 
zero and finite temperatures flSB&SBSHIIiEiiniE 
EQEIMESHIIll-At zero temperature a single mode 
of the system (the condensate) is macroscopically occupied 
and its evolution is well-described by the Gross-Pitaevskii 
equation. The classical field approximation is also suited to fi- 
nite temperature regimes where many modes of the system are 
highly occupied (also see j20ll ). The major advantage of clas- 
sical field approaches over other methods is that interactions 
between the modes can be treated non-perturbatively, making 
the formalism suitable for non-equilibrium studies. 

The purpose of this paper is to fully characterize the size 
and nature of the classical region for typical experimental pa- 
rameters. These properties of the classical region are not a 
priori obvious, especially for the case of an interacting gas. 
Better characterization of the classical region will be useful 
for: (a) informing better choice of parameters for classical 
field simulations of finite temperature gases; and (b) demon- 
strating the regimes of applicability of this theory. 



A. Historic Motivation: Blackbody Radiation 

We begin by reviewing a topic of modern physics that was 
responsible for the quantum revolution, yet also serves as 
motivation for the application of the classical field approxi- 
mation to quantum fields: the spectrum of radiation from a 
blackbody. Approximately a century ago the Rayleigh-Jeans 
theory of the light spectrum emitted from a blackbody was 
developed upon a classical field theory (i.e. Maxwell's elec- 
tromagnetic theory) in conjunction with statistical methods. 
The famous failure of this theory eventually led to the dis- 
covery that light is constituted of photons, and initiated the 
quantum revolution in physics. However, an observation of 
prime importance is that the Rayleigh-Jeans theory is actu- 
ally extremely successful at predicting the behaviour of the 
blackbody spectrum for the long wavelength modes. We il- 
lustrate this point in Fig. [Ha), where the quantum and classi- 



cal predictions for the spectrum are compared for a blackbody 
at a temperature of T = 2500 K. The most distinctive fea- 
ture of this graph is the rather alarming disagreement at small 
wavelengths, the so-called ultra-violet catastrophe. However 
at larger wavelengths, greater than 5/j.m say, the agreement 
between the quantum and classical predictions becomes in- 
creasingly good. For these long wavelength modes each pho- 
ton carries only a small amount of energy compared to k B T, 
so that on average these modes contain a large number of 
photons. This high occupancy is what we will refer to as 
the classical limit for a field. For these modes the discrete- 
ness of the energy each photon carries is masked by the large 
number of other photons in the same mode and the classi- 
cal theory provides a good description. We note that this is 
quite different from the high temperature classical limit for 
particles, when Bose/Fermi statistics can be approximated by 
Maxwell-Boltzmann statistics. The ultra-violet catastrophe 
occurs because the classical theory fails to describe the short 
wavelength modes correctly. These modes are sparingly oc- 
cupied by photons and require a proper quantum treatment. 
In Fig. 01b) the mean number of photons per mode is shown, 
and comparison with Fig. da) clearly reveals that the quan- 
tum and classical theories agree well where the mean number 
of photons is large. In general any electromagnetic mode can 
be made classical by making the temperature sufficiently large 
as photon number is not conserved. 



B. Classical Region of a Bose Gas 

Atoms are necessarily conserved and require a chemical po- 
tential (//) in their statistical description. By expanding the 
Planck and Bose-Einstein distributions for high occupation 
n(e) 3> 1 we find the following requirements on the single 
particle energies 



e < k B T 
e — fj, -C k B T 



Planck, (1) 
Bose — Einstein. (2) 



We note that in these limits the classical equipartition theo- 
rem holds l25ll . For the Bose gas at high temperatures (the 
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Figure 1 : (a) Quantum and classical theories for the energy spectrum 
of a blackbody at T = 2500 K. (b) The mean number of photons in 
each mode at T = 2500 K calculated using the Planck distribution. 



Maxwell-Boltzmann limit) the chemical potential is large and 
negative 



V (2irmk B T 
N 1 h? 



3/2 



(3) 



where V is the system volume and N is the number of atoms. 
In contrast to the photon case, this behaviour of /j, prevents the 
mean occupation of any individual mode from becoming large 
as T increases. 

Indeed, for a gas of atoms at temperatures above the mi- 
crokelvin regime, the average occupation of the system modes 
is much less than unity and the particle-like behaviour of the 
system dominates over the wave-like behaviour. However, at 
temperatures approaching the condensation temperature (T c ), 
the chemical potential approaches the ground state energy, and 
many modes of the system may become highly occupied. The 
nature of these highly occupied modes (outside of the con- 
densate itself) is not widely appreciated, and is the central 
topic we address in this paper. We refer to these modes as the 
constituting the classical region, for which the classical field 
approximation is appropriate: We can replace the quantum 
mode operators {a,j}, satisfying the commutation relations 



o», a 



Sij, with the c-number amplitudes {cj}. Modes 
outside the classical region are more sparsely occupied and 
are poorly described by the classical field approximation — 
these modes constitute the incoherent region. In this paper 
we characterize the classical region of a trapped Bose gas as 
a function of temperature and total particle number, and show 
that the classical field approximation should be widely appli- 
cable to current experiments. 

Schematically we show the classical and incoherent regions 
in Fig. [2] for the case of the harmonically trapped gas. The 
classical region is quantitatively defined as those single parti- 
cle modes with mean occupation greater than n cut , where n cut 
is the minimum occupation we require for a mode to be desig- 
nated as classical. As the single particle occupation monoton- 
ically decreases with increasing energy, the boundary of the 
classical region occurs at an energy we define to be the cut- 
off energy, e cu t. With reference to the electromagnetic case 
shown in Figs. Qla)-(b), we see that an adequate condition 
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Figure 2: Schematic diagram showing the classical and incoherent 
regions of the single particle spectrum for a harmonically trapped 
Bose gas. At the energy e cut the average occupation is n cut ~ 3 (see 
the text). 



for the classical field description to be valid is found by tak- 
ing n cu t ~ 3 (i.e. blackbody modes with a mean occupation 
of > 3 photons are equally well-described by the Planck and 
Rayleigh- Jeans distributions). We note that other authors have 
suggested that taking n cut ^ 5 to 10 particles may be more 
suitable @|. The precise value of n cut used will be unimpor- 
tant for the qualitative characteristics of the classical region 
we investigate here. 



C. Projected Gross-Pitaevskii Equation Formalism 

In this paper we are not primarily concerned with the details 
of applying the classical field technique to model Bose gases, 
but will briefly review the formalism to establish the context 
and importance of the classical region. The evolution equation 
for the classical modes of an interacting Bose field is give n by 
Projected Gross-Pitaevskii equation (PGPE) ifl l9i fTotllll : 



lh -m = 



where ^> = ^(x, t) is the classical matter wave field (i.e. de- 
scribes the atoms in the classical region), Vtrap(x) is the exter- 
nal trapping potential and Uq = AmaY? /m, with a the s-wave 
scattering length. The projector is defined as 



P{F(x)} ee 5>„(x) / d 3 xV*„(x')F(x'), 



(5) 



where cp n (x) are eigenstates of the single particle Hamiltonian 

H sp = — ^;V 2 + Vtrap( x ) an d the summation is restricted to 
modes in the classical region. The action of V in Eq. (0 is to 
project the arbitrary function F(x) into the classical region. 
A description of the incoherent particles will also be required 
for quantitative comparison with experiments and a formalism 
for including this into the PGPE theory has been presented in 
Ref. 

Previous work has shown that the PGPE will evolve ran- 
domised initial conditions to thermal equilibrium Jjl |H Oil 
HH |2j]]. The equilibrium state that is reached depends on 
three input parameters: (i) The number of atoms in the classi- 
cal region, which we refer to as iVbeiow; (ii) The total energy 
content of the classical region E; (iii) The cutoff energy e cut 
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determining the size of the classical region (i.e. the number of 
classical modes Mb e i ow HU). 

In dynamical simulations of Eq. (01, interactions cause the 
system to rapidly thermalize and subsequently (using the 
ergodic hypothesis) time-averaging can be used to extract 
equilibrium quantities. A major challenge for classical field 
approaches has been to determine thermodynamics quantities 
that are derivatives of entropy, such as temperature. This was 
recently overcome in Ref. QUI], using formalism initially de- 
veloped by Rugh ll22tl . A remaining issue with using the clas- 
sical field approach is that the thermodynamic parameters for 
the simulation (such as the temperature, and the total number 
of atoms in the system, i.e. including those above the cutoff) 
are only apparent a posteriori. The characterization of the 
classical region we give here should facilitate making a rea- 
sonable a priori estimate of parameters for PGPE simulations 
and should be of great practical benefit. 

n. FORMALISM 



Of interest for us is to identify the properties of the classical 
region, in particular the number of classical modes (Mb e i ow ) 
and their combined occupation (iVbciow)- According to Eq. 
(© the energy at which the mean occupation is equal to ?i cut 
is given by 









Hcut ) . 







Integrating the density of states (|6]l up to this energy we obtain 
the number of particle states, 

Mbeiow = / depHF(e), (10) 
Jo 

which we refer to as the number of classical modes. Finally, 
the combined occupation of the classical modes is given by 

A^bdow = / rfe"BE(e)pHF(e)- (H) 
Jo 



For the results presented in Figs. [3] and [4] we consider 
(bosonic) rubidium-87 in an isotropic harmonic trapping po- 
tential with single -particle energy spectrum e = hu>(n x +n y + 
n z + |), where uj = 2ir x 100 Hz is the harmonic trap fre- 
quency (chosen to be comparable with typical experiments). 
To account for interaction effects we use the semiclassical 
Hartree-Fock (HF) density of states 



Phf(c) 



d 3 x.d 3 p 

{2-Khf ( 



|- + Vtrap(x) + 2U n c (x) 



(6) 

where the condensate density is given by the Thomas-Fermi 
approximation 



n c (x) = 



I C^(MTF 



Krap(x)), M-Vt»p(x) >0, 

0, flTF - Vtrap(x) < 0, 

(7) 

and ptf is the Thomas-Fermi app roximation to the conden- 
sate eigenvalue (e.g. see 11231 12411 ). The Thomas-Fermi ap- 
proximation is only valid for large condensate number, how- 
ever, for small condensate numbers the density of states is 
largely unchanged from the ideal gas case, and so for our pur- 
poses this approximation is valid for all temperatures. 

To study the statistical properties of this system at finite 
temperature we work in the grand canonical ensemble ther- 
mally occupying the particle states described by (O using the 
Bose-Einstein distribution function 



"-be(c) 



1 



o0e 



lz-r 



(8) 



where /3 = 1/ksT is the inverse temperature, z = is the 
fugacity, and the excitation energy e is measured relative to 
the chemical potential (^tf)- To calculate the thermodynamic 
properties as a function of the total number of particles Nt it 
is necessary to find the fugacity subject to the requirement that 
the mean number of particles is equal to Nt- This parameter 
must be adjusted carefully as the condensate number, conden- 
sate density and excitation spectrum all change with z. 



III. RESULTS 
A. Ideal Gas 

Typical results for the mode occupations are shown in 
Fig. [3a) for a system of Nt = 2 x 10 6 atoms. Here we 
discuss the ideal gas results, and refrain from considering the 
interacting results until the next subsection. 

Far above the transition temperature, T c ks 566 nK, the oc- 
cupation of all modes is small compared to unity and there 
is no classical region. At T « 640 nK the occupation of 
several modes exceed n cut = 3 atoms and a classical re- 
gion develops in the system. At temperatures equal to, or less 
than T c , the occupation of the condensate mode dominates 
the system, however a large number of other modes have oc- 
cupations exceeding n cut . At temperatures near T c of order 
a thousand single particle states participate in the classical re- 
gion. This can be seen from Fig. [3{a) where the where the 
dashed line representing the occupation of the 1000 th single 
particle state is above the dotted line marking n cut for a range 
of temperatures encompassing T c . States of lower energy re- 
main highly occupied over a broader temperature range (e.g. 
the e = 8.5huj mode occupation, represented by the dash-dot 
line in Fig. [3{ a X exceeds n cut over a broader range than the 
e = n.hhio mode). Thus the number of modes in the clas- 
sical region Afbciow varies significantly with temperature, but 
is largest at T c where the non-condensate modes have their 
maximum occupation. 

The classical regions for two particular modes are shown 
in Fig. Ob) as a function of temperature and total number of 
particles. In region a the mean occupation of the e = 17.5fiw 
state exceeds n cut . In the broader region j3 the mean occu- 
pation of the e = 8.5hoj state exceeds n cut . The left hand 
edges of the regions are independent of Nt (i.e. the bound- 
ary of the region is vertical). This is because the thermal 
cloud is saturated for T < T c , and the mode occupation only 
depends on temperature. We can quantitatively predict the 
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Figure 3: (a) Mean occupation of the condensate and selected ex- 
cited single particle states for both ideal (black curves) and interact- 
ing (grey curves) systems of 2 x 10 6 rubidium-87 atoms in a 3D 
100 Hz isotropic trap. Condensate: dash-dot. M = 120 th single 
particle level: solid. M = 1000 th single particle level: dashed. 
The horizontal light-grey line indicates n cu t = 3 particles per mode. 
The vertical light-grey line marks T c for this system, (b) Phase di- 
agram of the classical region. The dark shaded region labelled a is 
where the M = 1000 th single particle level has a mean occupation 
of greater than 3 particles. The lighter shaded region labelled j3 is 
where the M = 120 th single particle level has a mean occupation 
of greater than 3 particles. The boundaries of these regions for the 
non-interacting case are shown as grey lines and T c is shown as with 
black dots. 



left-hand boundary by noting that below transition tempera- 
ture the chemical potential is well approximated by taking it 
equal to the ground state energy, i.e. fx « e , inverting the 
Bose-Einstein distribution function we obtain 



Tlh — 



e-e q 



k B ln(l + 1/ricut) ' 



(12) 
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independent of Nt and in agreement with the left-hand 
boundaries in Fig.[3jb). 



Figure 4: The classical region for rubidium-87 atoms in a 3D 
isotropic harmonic trap with u — 2n x 100 Hz. (a) The total num- 
ber of atoms in modes with occupation greater than 3 (iVbeiow) for 
systems with a total of 5 x 10 5 (solid line), 1 x 10 6 (dash-dot line), 
and 5 x 10 6 (dashed line) atoms. The dotted curves show the number 
of atoms constituting iVbeiow that are not in the condensate, (b) The 
number of modes in the classical region Afbciow for the systems in 
(a). For all curves in this plot black curves are for the non-interacting 
case and grey curves include the effects of interactions. For reference 
T c is indicated using vertical dotted lines in each plot. 



The right hand boundary occurs above T c where [i is differ- 
ent from eq, and is seen to be dependent on Nt- We note 
that for large Nt the classical region extends to tempera- 
tures considerably above T c , indeed for the case considered 
in Fig. [2b) the classical region begins ~ 100 nK above T c for 
N T > 3 x 10 6 particles. 

In Fig.[4]we characterise the total number of atoms (A^iow) 
and modes (Afboiow) participating in the classical region as a 
function of temperature for several values of Nt- For T < T c , 
-Wbeiow is dominated by the condensate mode and scales as 
-^beiow ~ {1 — {T/T c ) 3 }, as can be seen for the three cases 
of Nt considered in Fig.[4ta). The number of classical region 
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particles not in the condensate is shown as the dotted curve in 
Fig. Ufa). We notice that up to the transition temperature this 
curve is independent of Nt, due to the saturated nature of the 
thermal cloud. Though the population of these classical non- 
condensate particles is much less than the condensate, for the 
cases considered in Fig. Hta) we see that they may constitute 
up to ~ 10 5 particles. Additionally, above T c there is no con- 
densate and the classical non-condensate particles dominate 
the classical region. 

In Fig.Htb) we consider the behaviour of the cutoff energy. 
For T < T c the cutoff energy increases linearly with tem- 
perature, and is independent of Nt- As discussed above, for 
T < T c the chemical potential is well-approximated as the 
ground state energy fj, — ► eo, and we can find an explicit ex- 
pression for the cutoff energy in terms of the temperature 



e cu t = k B T\n ( 1 



1 

n cu t 



eo, 



(T < T c ), (13) 



(or equivalently, A/bciow = e cut/6(^0 > where u> is the geo- 
metric mean of the harmonic trap frequencies). Reaching its 
maximum value at T c , the cutoff energy decreases for T > T c 
as the classical region shrinks due to the rapid decrease in 



B. Hartree-Fock Results 

The results for the interacting gas, calculated using the den- 
sity of states given in Eq. ©, are also shown in Figs. [3] and |U 
There are many qualitative differences introduced in the inter- 
acting case. First, for the interacting system the size and num- 
ber of particles in the classical regions tends to be larger than 
the ideal case. Second, the number of modes of the classical 
region (Mbeiow) tends to peak at temperatures below T c for 
the interacting case, whereas it peaks at T c for the ideal case. 
These effects are related and can be understood by examining 
the density of states (O. Due to the mean field repulsion of 
the condensate, the lowest energy excitations do not occupy 
the phase space coordinates near the origin, but for position 
values near the surface of the condensate. Because of the in- 
creased amount of phase space available at the surface of the 
condensate more low energy states are available. As the tem- 
perature increases the condensate fraction and radius both de- 
crease, reducing this phase space enhancement. The observed 
maximum of Mbeiow in Fig. Hfb) occurs due to the competi- 
tion between increasing temperature (which tends to increase 
the number of accessible states) and reducing condensate ra- 
dius (which reduces the number of low energy modes). 



IV. COMPARISON TO PGPE 

It is of interest to assess how well the ideal gas and Hartree- 
Fock results for the classical region describe the classical field 
case. Because the PGPE approach includes interactions be- 
tween the excitations we expect quantitative differences with 
the Hartree-Fock results to arise. Our procedure for calcu- 
lating equilibrium properties for the trapped Bose gas is dis- 




Figure 5: Comparison of classical region estimates using the meth- 
ods described in this paper to the numerical PGPE results as a 
function of temperature, (a) Number of classical region modes for 
n cu t = 2.8: HF (solid), ideal gas (dashed), PGPE solution (+), in- 
dicating estimated error bars vertically, (b) Classical region frac- 
tional population: HF (dotted), ideal gas (solid), PGPE solution 
(square). Condensate fraction: HF (dashed), ideal gas (dash-dot), 
PGPE solution (circle). Harmonic trap with oscillation frequencies 
{fx, fy,.fz} = {1, 1, \/8} x40 Hz containing approximately 3 x 10 5 
87 Rb atoms. 



cussed in Ref. @]. Briefly, the PGPE system is a microcanon- 
ical system with external constants of motion (i.e. thermo- 
dynamic constraints) of N^eiow an d total energy E. Addi- 
tionally, to avoid an ultraviolet divergence a restriction in the 
single particle modes must be made by selecting Mbeiow- Af- 
ter equilibrium is found and the properties of the classical re- 
gion are determined, such as mean density, temperature and 
chemical potential, the incoherent region can be found using 
a semiclassical Hartree-Fock analysis. 

The results of this study are shown in Fig. [5] These confirm 
that the HF predictions are qualitatively in much better agree- 
ment with the PGPE theory than the ideal case. We note the 
following points: 

1. Temperature in Fig. [5] has been scaled in terms of the 
ideal transition temperature. The PGPE exhibits the 
condensation transition at T/T c < 1 due to the down- 
ward shifts in critical temperature arising from finite 
size and interaction effects (neither of these effects are 
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in the ideal or HF results). 

2. The PGPE results suggest that number of classical 
modes, Mb e iow> does not increase up to the critical point 
(as predicted by the ideal result), and are consistent with 
the HF prediction of a maximum value occurring for 

T < T c . 

3. The number of classical region modes is considerably 
larger than predicted by the HF method. However, as 
the total number of modes Mbciow is proportional to 
ejjjjt for the ideal gas then the increase in e cut is not as 
significant. 

V. CONCLUSIONS AND OUTLOOK 

We have characterized the classical region, containing 
modes with occupations higher than n cut ~ 3, as a function 
of the total number of particles and temperature for a harmon- 
ically trapped Bose gas. To do this we have examined the 
number of modes and number of particles in the classical re- 
gion, and the behaviour of the energy cutoff. 

For the ideal case we have shown that for T > T c the cut- 
off energy is dependent on the total number of atoms in the 



system and the temperature, whereas for T < T c the cutoff 
energy increases linearly with temperature and is independent 
ofN T . 

We have studied interaction effects using a Hartree-Fock 
approach which we have compared with full PGPE simula- 
tions. Our results show that interactions have a considerable 
effect on the classical region, and appear to make the size of 
the classical region much greater than predicted by an ideal 
gas analysis. 

In general we have demonstrated that for systems with of 
order a million particles the classical region may consist of 
thousands of single particles states and extend to temperatures 
of order 100 nK above T c . Overall these results demonstrate 
that the PGPE technique has a wide range of applicability to 
trapped Bose gases and should facilitate better parameter es- 
timation necessary for effectively applying the PGPE formal- 
ism to model experiments. 
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